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I.  Introduction 

Intuitively,  if  one  has  two  sets  of  environmental  data  which 
are  "close,"  and  one  uses  them  as  inputs  to  a  transmission  loss 
model,  then  one  would  expect  the  resulting  transmission  loss  curves 
to  be  "close."  Furthermore,  if  two  sets  of  environmental  data  have 
little  or  nothing  in  common,  one  is  not  surprised  if  their  resulting 
transmission  loss  curves  vary  greatly.  From  this  rather  modest 
hypothesis,  one  can  develop  some  sort  of  mental  "scale"  which 
qualitatively  compares  transmission  loss  curves.  Such  a  scale  could 
have  values  such  as  "seems  close,"  "looks  entirely  different,"  and 
any  number  of  other  shades  of  grading.  Yet,  such  a  concept,  often 
based  on  visual  inspection,  is  still  very  subjective. 

If  one  is  going  to  seriously  consider  the  U.S.  Navy's  ability 
to  predict  transmission  loss,  then  it  is  inadequate  to  merely  address 
the  models:  transmission  loss  models  currently  exist  which  are  quite 
sophisticated.  One  must  also  consider  the  variability  (or,  in  some 
cases,  the  accuracy)  of  the  environmental  inputs.  The  link  between 
environmental  changes  and  transmission  loss  changes  is  the  model 
itself.  If  one  is  to  analytically  approach  questions  dealing  with 
similarities  or  differences  in  transmission  loss  curves  produced 
from  different  input,  it  is  necessary  to  establish  quantitative 
measures  on  the  transmission  loss  curves:  "close"  is  not  good 
enough.  The  purpose  of  this  paper  is  to  develop  a  series  of  measures 
which  may  be  applied  to  transmission  loss  curves  to  quantitatively 
describe  them  and  afford  a  quantitative  means  of  comparison. 

It  should  be  noted  at  the  outset,  however,  that  the  intuitive 
statement  opening  this  report  is  not  always  true.  A  small  change  in 
the  environment  can  dramatically  alter  the  transmission  loss  at  a 
point  (  cf.  Ch.  II)  by  shifting  a  caustic.  And,  if  the  sound  speed 
is  altered  so  that  the  change  in  the  profile  c(z)  is  small,  but 
the  change  in  the  derivative  with  respect  to  depth,  3c/3z,  is  large, 


then  the  transmission  loss  field  can  be  very  different.  But  "in 
general,"  one  will  observe  that  a  small  change  in  the  input  usually 
results  in  a  small  change  in  the  transmission  loss  curve. 

In  regard  to  our  approach,  one  would  think  that  a  metric1  to 
measure  "closeness"  of  two  transmission  loss  curves  could  be  devised, 
i.e.,  given  any  two  transmission  loss  curves,  one  could  plug  them 
into  a  formula  and  come  up  with  a  single  number.  This  is  indeed 
true,  but  what  meaning  would  it  have?  Consider  a  transmission  loss 
curve  which  is  basically  10  log  r  loss,  with  prominent  convergence 
zones.  Does  one  say  that  this  curve  is  "closer"  to  an  18  log  r 
loss  curve  with  the  same  prominent  convergence  zones  than  it  is  to  a 
monotonic  10  log  r  loss  curve  with  no  convergence  zones?  One  is 
essentially  comparing  apples  and  oranges,  and  it  is  felt  that  such 
metrics  tend  to  confuse  rather  than  enlighten.  Thus,  we  have 
developed  a  series  of  measures  describing  various  characteristics 
of  the  curves. 

II.  Statement  of  Problem 
2 

Given  two  transmission  loss  curves  (i.e.,  TL  in  dB  re  some 
standard  versus  range  in  some  units) ,  define  some  quantitative 
measures  to  indicate  the  similarity  or  difference  of  the  curves. 

These  measures  must  be  such  that  a  computer  may  obtain  them,  and  do 

,  3 

so  m  a  reasonably  time-efficient  manner. 

The  approach  of  looking  at  transmission  loss  verses  range  means 
that  curves  of  the  same  environmental  parameters,  but  differing 
receiver  depths,  may  be  compared.  But,  it  precludes  the  direct 
comparison  of  an  entire  transmission  loss  field  (i.e.,  a  grid  of 
TL  at  various  ranges  and  depths)  with  another  transmission  loss 
field. * 

Further  difficulties  arise  through  the  mode  of  describing  the 
transmission  loss  curves.  For  although  the  curves  theoretically 
extend  over  all  positive  reals,  in  actuality  one  does  not  measure 
transmission  loss  to  an  infinite  distance  nor  to  an  infinitesimal 

^Since  footnotes  have  an  occasional  but  unfortunate  habit  of  being 
lengthy,  they  have  all  been  relegated  to  a  section  in  the  rear, 
following  the  appendices.  Note  1  will  be  found  there. 


/.'V- 
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distance.  Thus,  the  curve  whose  domain  is  (0,  °°)  ,  is  actually 
specified  on  a  closed  subset  of  this  open  interval.  Furthermore, 
whether  the  transmission  loss  is  measured  or  calculated,  it  is 
represented  not  by  a  curve,  but  by  a  discrete  set  of  points. 

When  one  compares  two  transmission  loss  "curves"  in  a  computer,  one 
is  actually  dealing  with  two  sets  of  discrete  points,  and  their 
grids  may  be  completely  different.  This  implies  that  simple 
expressions  such  as  the  difference  of  two  transmission  loss  curves 


T1 “  T2  (r) 


may  not  be  directly  computable  from  the  data.  It  will  always  be 
assumed,  however,  that  we  are  working  with  transmission  loss  curves, 
not  points,  as  curves  may  be  constructed  by  interpolation  between 
the  points.5 

Now  that  the  problem  has  been  stated,  along  with  a  few  rele¬ 
vant  observations,  one  may  briefly  consider  the  classical  approach. 
Let  T-^r)  and  T2  (r)  be  two  transmission  loss  functions,  and  assume 
that  the  closed  interval  [a,  b]  is  within  the  domain  of  each  func¬ 
tion.®.  For  any  pil,  we  may  consider  the  L  norm  of  the  differences, 

hr 

i.e. 


T1  "  T2  I'  p  =  lP<t1-t2>  = 


.  ,b 

5=5-  J  |T1(r)-T2(r)  |  pdr 


1/P 


The  three  common  norms  are 

P 

1 1 T1  -  T2  ||  ^  =  Average  difference  between  the  curves 

HT1  ”  t2  II  2  =  Root-Mean“Sc3uare  deviation  of  curves 

||  T^  -  T2  II  oo  =  Maximum  difference  between 

These  measures,  for  p>l,  are  quite  handy  and  are  often  used7,  but 

they  do  not  provide  an  adequate  measure  for  the  situation  where 
the  transmission  loss  curves  differ  by  the  horizontal  shift  of  a 


3 


convergence  zone.  Suppose  the  curve  has  a  caustic  followed  by 
a  shadow  zone,  with  a  drop  of  say  30  dB.  If  a  slightly  different 
environment  is  used,  a  new  curve,  T-  may  be  produced  which  is 

^  O 

similar,  buth  with  a  slight  horizontal  shift  (Fig.  1).  Then, 
although  the  norm  is  still  small,  both  and  norms  of  the 
difference  will  be  large.  Since  shifting  of  convergence  zones  is 


quite  common  when  a  Sound  Speed  Profile  is  altered,  it  is  seen 


that  point  criteria7  (such  as  Lp  norms)  are  not  adequate  for  our 


needs. 


III.  Quantitative  Characterization  of  Curves 


A  typical  transmission  loss  curve  (Fig.  2)  may  be  described 
(qualitatively  or  quantitatively)  in  terms  of  its  typical  character 
istics : 


i) 


Monotonic  Element:  If  one  were  to  "smooth  out"  all  the 
convergence  zones,  so  that  the  loss  was  principally 
spreading  and  attenuation,  then  what  loss  curve  does 
one  obtain? 


ii) 


Existence  of  Convergence  Zones:  Most  long-range  propa¬ 
gation  loss  functions  show  zones,  but  there  are  excep¬ 
tions.  If  there  are  convergence  zones,  then  one  may 
consider  the  next  two  characteristics. 


iii)  Convergence  Zone  Placement:  This  concerns  the  distance 
from  the  source  to  the  first  convergence  zone,  and  the 
distances  between  the  successive  convergence  zones. 

The  distance  between  successive  zones  need  not  be  con¬ 
stant,  but  in  most  cases,  it  will  not  vary  radically. 

iv)  Convergence  Zone  Shape:  This  characteristic  includes 
properites  such  as  the  geometric  shape  of  the  zone  (on 
both  sides  of  the  caustic) ,  the  spread  of  the  successive 
convergence  zone  widths,  and  the  decay  (or  growth)  of  the 
successive  convergence  zones.  Properties  of  the  zones 


Figure  2.  Typical  TL  curve  showing  descriptive  elements 


such  as  the  difference  in  intensity  between  a  conver¬ 
gence  zone  and  the  adjoining  shadow  zone  might  also  be 
considered  as  characteristics  of  shape,  but  they  will 
be  dealt  with  when  considering  deviations  from  the 
monotonic  component. 

Each  of  these  may  be  quantitatively  described  by  a  variety  of 
methods,  depending  upon  the  amount  of  detail  with  which  one  wishes 
to  work.  An  inadequate  level  of  detail  results  in  phenomena  being 
lost  through  simplification,  while  an  excessive  level  of  detail 
results  in  phenomena  being  buried  in  myriads  of  numbers.  Neither 
extreme  is  acceptable.  The  approach  given  below  is  based  on  analytic 
considerations  in  these  four  categories  of  characteristics. 


IV.  Definitions 

In  this  section,  formal  definitions  are  given  for  a  number 
of  variables  and  functions  which  will  be  used  throughout  the  paper. 
For  quick  reference,  a  list  of  notation  is  also  provided  following 
the  appendices. 

Let  Tk(r),  k  =  1 ,  2  be  transmission  loss  curves,  which  are  given 
out  to  range  r^.  One  could  consider  that  the  domain  of  Tk  was  the 
half-open  interval  (0,  r^],  but  transmission  loss  data  is  not 
collected  down  to  zero  range,  nor  are  long-range  models  always 
accurate  at  absurdly  short  ranges.  Thus,  we  choose  a  fi>0  and  con¬ 
sider  the  transmission  loss  functions  to  begin  at  range  <$.  Let 
r^2  be  the  minimum  of  r^  and  r2  and  we  will  always  assume 


Thus,  the  domain  of  Tk  will  be  the  closed  interval  [6,  r^] .  Let 
£k(r)  be  a  monotonic  loss  function  corresponding  to  Tk(r),  so  that 


its  domain  is  also  [6,  r^] . 


Let  £^2  (r)  be  a  monotonic  loss  func¬ 


tion  corresponding  to  T^  and  T where  domain  will  be  [6,  r123* 


(The  derivation  of 
And,  let  Pk(r) ,  Pkk(r) 


and  f°un<3  in  the  next  section.) 

be  the  perturbation  components,  defined  by 


Pk(r)  =  Tk(r)  "  £k(r) 
Pkk(r)  =  Tk(r)  "  ^12  (r) 


As  discussed  in  Chapter  X,  the  perturbation  component  p  (r)  will 

not  normally  be  used;  we  include  -  here  for  convenience  of  defini¬ 
tion.  It  should  also  be  noted  t  t  when  functions  of  range  (r)  are 
explicitly  stated,  it  is  assumed  aat  r  is  in  nautical  miles. 

Now,  we  turn  to  the  develop;  '  .  of  a  weighting  function  for 
integration.  Let  f  (r)  be  a  weighting  function  defined  on  (0,  °°) 
which  will  be  used  to  ensure  that  the  analysis  is  not  strongly 
dependent  upon  r^.  Ideally,  f  should  be  more-or-less  level  across 
the  range  of  direct  propagation,  first  convergence  zones,  first 
shadow  zone,  and  second  convergence  zone.  It  should  then  drop  off 
and  approach  zero  asymptotically.  Three  typical  forms  for  f  are 
shown  in  Fig.  3. 

The  linear  form 


1 

for 

r  <  70 

110-r 

40 

for 

70  <  r  <  100 

1/4 

for 

100  -  r 

is  very  convenient  for  computer  calculation,  but  is  a  somewhat  crude 
approximation,  particularly  in  regard  to  the  asymptote.  .  The  rational 
form 


f  (r)  =  l/(l+x4) 


where 


b)  Rational  Approximation 


c) 

Figure 


70  nm  90  nm  110  run 


Ideal  form 

3.  Three  forms  for  the  weighting  function. 


is  vastly  superior  in  this  regard.  A  more  precise  representation 
of  the  weighting  function  (Fig.  3c)  may  be  obtained  using  hyperbolic 
trigonometric  functions,  but  the  complexity  of  the  computer  calcula¬ 
tions  exceeds  the  marginal  benefit  gained  over  the  rational  approxi¬ 
mation. 


It  was  stated  above  that  the  shape  of  f  was  chosen  so  that  the 
analysis  would  not  be  strongly  dependent  upon  the  values  of  the 
two  r^.  This  statement  will  now  be  explained.  The  analysis  used 
in  this  paper  is  based  on  various  integrals  of  the  transmission 
loss  functions.  These  integrals  are  usually  algebraic  combinations 
of  the  functions,  hence,  it  suffices  to  consider  the  most  elementary, 
i.e.,  a  simple  average. 


i 


T1  (r) dr 


y(r1)  = 


It  is  apparent  from  Fig.  4  that  this  average  is  highly  dependent 
upon  the  argument,  for  y (v^  differs  considerably  from  y(v2).  But, 
the  analysis  of  propagation  loss  curves  should  not  be  unduely 
influenced  by  the  amount  of  the  curve  (i.e.,  the  range)  which  is 
known. Therefore,  integrals  such  as  the  simple  average  given 
above  are  not  a  good  foundation  upon  which  to  build  the  theory. 

To  avoid  such  a  situation,  a  weighting  function  is  introduced  which 
goes  to  zero  asymptotically  as  the  range  increases.  Since  we  are 
dealing  with  long-range  propagation  modeling,  one  may  expect 
that  most  propagation  loss  curves  will  be  known  to  ranges  over  70  nm, 
and  that  a  range  of  110  nm  is  often  attained.  Thus,  in  this  region, 
the  weighting  function  should  decrease  substantially.  One  usually 
encounters  the  first  and  second  convergence  zones  and  the  first 
shadow  zone  within  70  nm  or  so  of  the  source,  so  that  these  primary 
characteristics  will  be  weighted  heavily.  Conversely,  many  model 
runs  are  made  to  several  hundred  miles.  But,  if  one  is  comparing 
two  propagation  curves,  it  is  not  clear  that  one  could  justifiably 
compare  the  finer  properties  of  the  curves  at  ranges  of  several 
hundred  miles.  Thus,  the  low  weight  given  to  these  ranges  is 
acceptable. ^ 


V-** 
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V.  The  Monotonic  Component 

In  this  section,  we  investigate  the  spreading  and  attenuation 
properties  of  the  loss  curve.  Let  us  assume  that  these  two  phenomena 
may  be  described  by  a  function  of  the  form 

(r)  =  b,  +  b.  log  r  +  ar 


where  bk  b^  are  constants  to  be  determined  and  a  is  the  attenuation 

coefficient. ^ 3  For  any  function  £  of  suitable  domain,  define  the 
functional 


0kU)  = 


t  K  2 

J  (Tk(r)-£(r)rf  (r)dr 

/  kf (r) dr 
6 


Then,  we  choose  bk,  bk  (k=l,  2)  such  that  £^  minimizes  0^  and  £2 

minimizes  02.  The  function  £k  represents  the  best  weighted  least 

squares  fit  from  this  class  of  monotonic  functions  to  the  curve  Tk> 

The  function  £k  will  be  used  as  the  descriptor  of  the  monotonic 

component  of  Tk.  In  Chapter  X,  which  deals  with  the  comparison  of 

two  propagation  loss  curves,  it  will  be  necessary  to  use  another 

curve  which  represents  the  spreading  and  attenuation  of  both  loss 

curves.  This  common  monotonic  component  will  be  a  function  £. » 

14  . 

of  the  above  class,  and  it  will  be  the  unique  function  which 
minimizes  the  functional 


012U)  =  01U)  +  02U) 


The  three  values 


Ql  (£l)  ,  ®2  U2)  ,  0i2  ^12^ 

are  statistics  which  should  be  computed  when  comparing  two  curves, 
An  efficient  method  for  obtaining  these  quantities  is  given  in 


*  wrv  ,*v 


Appendix  A.  These  are  measures  of  the  deviation  of  the  curves 
from  a  monotonic  loss  function. 

Another  quantity  of  interest,  to  measure  sudden  rise  or  drop 
in  the  transmission  loss  curve,  will  be  defined.  If  one  wished  to 
measure  such  a  characteristic,  then  the  difference  would  be 

lTk(r)  "  Tk(r+0)  I 

where  o  was  some  small  amount  (e.g.,  several  miles).  Thus,  one 
might  consider  a  measure  such  as 

max{  |  T^  (r) -T^  (r+a)  |  :  6  £  r  £  rk~a,  0  <  a  r  p} 

where  p  is  a  constant  (e.g.,  10  nm) .  Unfortunately,  while  this  is 
a  very  good  measure,  it  has  the  one  drawback  of  being,  in  certain 
cases,  strongly  dependent  upon  r^.  Thus,  one  chooses  a  less  mean¬ 
ingful  but  more  stable  measure 

Xk  =  max{  |Tk(r)-Tk(r+a)-bk  log  (1+j-)  -ao  | yjf  (r )  f  (r+o)  : 

6  <  r  *  r k  -  a,  0*o£p} 

Let  ak  be  the  value  of  o  which  maximizes  the  functional. 

VI.  Convergence  Zones:  To  be  or  not  to  be 

The  decision  as  to  whether  or  not  a  transmission  loss  curve 

has  convergence  zones  would  appear  to  be  relatively  straightforward: 

an  acoustician  can  make  such  a  determination  by  visual  inspection 

(cf .  Fig.  5) .  But,  such  a  determination  may  be  difficult  when  one 

has  points  arising  from  either  measured  data  or  a  coherent  model 

which  exhibit  a  significant  scatter,  so  that  the  presence  or  absence 

15 

of  slight  convergence  zones  is  hidden  by  the  "noise"  (Fig.  6). 


Figure  5.  Obvious  presence  and  absense  of  convergence 


One  needs  a  method  of  discrimination  which  the  computer  can 
perform,  and  a  three-step  procedure  is  used.  To  determine  if  (r) 
has  convergence  zones,  the  first  step  is  to  examine  Q^(£^),  which  is 

a  measure  of  the  deviation  of  T^(r)  from  (r)  .  Such  deviation  may 

be  due  to  either  data  scatter  (cf.  Fig.  6)  or  to  convergence  zones. 

Let  ic  be  a  threshold  such  that  if 

©i ( ^! )  <  *2 

then  one  immediately  concludes  that  convergence  zones  are  not 
present.  Let  <  be  a  threshold  such  that  if 

k2  <  01(£1)  <  k2 

then  it  is  very  probable  that  convergence  zones  are  not  present, 

but  the  second  step  should  be  used  for  this  determination.  Reasonable 

values  for  <  and  k  might  be  2  dB  and  2.5  dB,  respectively. 

The  second  step  is  examination  of  xk  to  try  and  determine  whether 
a  small  amount  of  deviation  is  due  to  data  scatter  or  convergence 
zones.  This  is  only  done  when  0^(JI^)  lies  between  the  two  thresholds. 
If 

>  k2 

then  this  step  should  be  ignored  and  step  three  used  for  determina¬ 
tion.  If  T^(r)  was  identical  to  its  monotonic  component,  then 

would  be  zero.  Thus,  one  may  choose  a  threshold  0  and  state  that  if 

X2!  <  0*0!  Ux) 

then  the  variation  of  T^  from  Jl^  is  a  slowly-varying  function;  too 

slowly  varying  to  be  convergenze  zones. ^  This  threshold  is  very 
difficult  to  estimate;  a  starting  value  might  be  1/3.  And,  a  thresh¬ 
old  d<p  is  used  to  test  for  short-range  fluctuations.  If 


then  is  assumed  that  we  are  dealing  with  short-range  fluctuations 
rather  than  convergence  zones.  Ideally,  we  should  have 

o  <<  ~p  . 

But,  the  transmission  loss  curve  is  entered  at  discrete  points, 

and  from  the  definition  of  x^/  it  is  obvious  that  will  be  either 
p  or  the  distance  between  two  of  these  discrete  points.  For  short- 
range  fluctuations,  will  not  be  p  (cf.  Fig.  21a),  so  o1  will  be 

a  distance  between  some  pair  of  input  points.  Thus,  the  test  with 
o  is  nontrivial  only  if  5  is  greater  than  the  distance  between 
adjacent  points  of  the  tabulated  transmission  loss  data.  If  trans¬ 
mission  loss  data  was  given  every  10  km (=6. 2  nm)  and  p  is  10  miles, 
then  it  is  obvious  that  we  cannot  have  both 

6.2  i  5  <  10 

and 

5  <<  10  . 

Therefore,  this  criteria  may  be  difficult  to  successfully  implement 
if  the  grid  spacing  of  the  tabulated  transmission  loss  is  large.  If 
is  greater  than  that  threshold,  then  one  would  proceed  to  step 

three. 

The  third  step  is  to  determine  whether  the  suspected  conver¬ 
gence  zones  really  are  significant  zones.  This  step,  which  is  the 
final  one  of  the  convergence  zone  determination,  is  done  in  a  con¬ 
structive  manner;  i.e.,  we  find  the  first  convergent  zone.  If  the 
algorithm  to  find  the  first  zone  fails,  one  assumes  that  no  con¬ 
vergence  zones  exist.  This  algorithm  is  described  in  the  next  sec¬ 
tion,  Convergent  Zone  Placement.  If  no  such  zone  is  detected,  then 
the  characterization  of  the  propagation  loss  curve  is  limited  to  its 


monotonic  component,  the  mean  square  deviation  from  that  component, 
and  the  maximum  variation  over  an  interval.17 


VII.  Convergence  Zone  Placement 

In  this  section,  we  restrict  attention  to  the  problem  of 
determining  the  location  of  each  convergence  zone  (the  peak,  not 
the  full  zone).  From  this,  we  will  develop  measures  which  can  be 
used  to  characterize  convergence  zone  placement. 

An  obvious  way  to  obtain  the  distance  to  the  initial  zone, 

and  the  spacing  of  successive  convergence  zones,  is  to  take  the 

(local)  maxima  (Fig.  7).  In  actual  use,  this  may  be  done  easily 

by  eyeball,  but  not  by  the  computer,  for  many  transmission  loss 
1 8 

curves  are  very  jagged  (Fig.  8). 

Merely  identifying  the  maximum  value  within  a  range  of  one 
convergence  zone  is  possible,  but  this  will  lead  to  erroneous  identi¬ 
fication  in  certain  instances  (e.g..  Fig  9).  A  straightforward 
solution  to  this  is  to  use  a  running  weighted  average  of  the  trans¬ 
mission  loss  curve,  not  the  curve  itself.  Then,  the  local  maxima 

of  the  weighted  sliding  average  will  correspond  to  the  convergence 
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zone  peaks.  To  obtain  such  a  running  average,  choose  a  function 
coR(r)  such  as  that  in  Fig.  10  as  the  weighting  kernel.  Appendix  T 

discusses  typical  forms  for  wR(r)  •  Using  this  function,  define 

*k(R)  *  Pk(r)uR(r)f  (r)dr 
as  the  averaged  curve. 

Since  the  problem  is  to  identify  the  first  peak,  a  range  Rq 

which  is  past  the  (expected)  position  of  the  peak  is  chosen,  and 
$k(R)  is  examined  on  the  interval  [6,  RQ] .  This  function  should 
look  more-or-less  like  that  of  Fig.  11.  The  desired  solution  is  a 
local  maxima, ^  i.e.,  a  point  where  ^'^(R)  changes  sign  and  $"^(R) 

is  negative.  This  value  of  R  will  be  the  location  of  the  initial 
peak;  call  it  R^.  When  it  is  important  to  designate  it  as  the 
initial  peak  of  Tu (r) ,  it  will  be  denoted  R^  One  could  continue 


One  could  continue 


Figure  9.  situation  in  which  actual 
maximum  of  TL  curve  is  not 
best  estimate  of  convergence 
zone. 


w  (r) 
R 


Figure  10.  Weighting  function 


but  such 


this  method  to  locate  the  successive  peaks  ,  R^,  etc., 

a  method  is  fraught  with  numeric  instabilities  which  may  easily 
result  in  false  "convergence  zones." 

It  should  be  pointed  out  that  one  may  have  a  peak  where  the 

2  * 

transmission  loss  is  still  monotonic,  as  in  Fig.  12. 

From  Appendix  B  it  is  seen  that  the  spacing  between  convergence 
zones  is  not  (in  general)  constant,  but  the  deviation  from  a  constant 
is  rather  small.  Therefore,  rather  than  compute  the  intervals 
between  the- successive  convergence  zones,  the  average  interval 
between  convergence  zones  will  be  estimated  and  that  will  be  the 
descriptor. 

If  one  writes 


A  .  R  =  A  +  l  . 

3  D 

where  the  are  small  and  "average  out"  to  zero,  then  A  will  the 
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the  desired  spacing  estimate.  To  this  end,  define  the  function 


°A(r)  u)R1+jX(r) 


j=0 


for  X>0.  Then,  the  value  Xk  which  is  a  local  maximum  for  the 
functional 

=fZk  Pk(r)Qx  (r)f(r)dr//rkQx  (r)f(r)dr 


will  be  the  estimate  of  A  for  the  transmission  loss  curve  Tk(r) 
Since 
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6  *  Ri  s  Ro  s  rk 
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the  denominator  of  this  function  cannot  be  zero.  '  Define 

nk ,«rl  rk  pk  <r>  %  <r> £  (r>  ^ ar//{rk  “Rl  V)  f  (r)dr 


13 


lim  ft(X,  ) 

xk.o 

2  6 

exists  and  except  for  some  pathological  cases,  it  is  less  than 

27 

oo.  So,  a  true  maximum  will  generally  exist,  and  Fig.  13  illus- 

2  8 

trates  a  common  case  for  which  this  is  true.  But,  Fig  13  also 
illustrates  a  case  where  the  true  maximum  is  not  the  desired  maximum. 
Figure  14  shows  another  case  where  the  true  maximum  is  not  the 
desired  maximum.^® 

It  is  clear  from  these  examples  that  when  Xk  is  chosen,  it 
should  not  be  a  global  maximum,  but  rather  a  local  maximum.  Thus, 
an  interval  must  be  chosen  so  as  to  span  only  one  convergence  zone 

,  3i 

at  a  time,  and  the  function  is  examined  on  such  an  interval. 

In  summary,  this  chapter  has  given  the  techniques  to  determine 
the  descriptors  k  and  X^  of  the  convergence  zone  placement  of 

Vr)- 


VIII.  Convergence  Zone  Shape 

In  this  section,  a  quantitative  measure  is  associated  with  the 
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shape  of  a  convergence  zone.  There  are  a  number  of  measures  one 
could  devise;  those  chosen  are  easily  implemen table.  It  is  antici¬ 
pated  that  the  other  measures  (e.g.,  monotonic  component,  deviation 
from  monotonic  component,  convergence  zone  spacing)  will  prove  more 
significant  in  practice. 


One  could  attempt  to  describe  the  shape  of  each  convergence 
zone,  or  describe  the  overall,  or  "average"  shape  of  the  convergence 
zones.  Or,  one  may  merely  describe  the  shape  of  a  selected  set  of 
zones  {e.g.,  the  first  two).  It  is  this  last  approach  that  will  be 
taken,  although  the  technique  used  is  basic  to  all  three  approaches. 

Let  Rj  be  the  (averaged)  distance  to  the  j-th  convergence  zone 
from  the  source,  i.e., 

Rj  =  rx  +  (j-mk 

and  consider  the  shape  of  the  j-th  peak.  The  shape  will  be  modeled 
by  the  function 

K  +  d2x  +  jtIt 
Vr>  =  d  ° 

K  +  a5x  +  ^S~ 

o 

where  d^  are  constants  and 

x  =  -r  -  Rj  . 

The  determination  of  these  coefficients  is  found  in  Appendix  A. 
Although  there  appear  to  be  seven  coefficients,  one  of  them  is 
chosen  beforehand,  and  the  condition  is  imposed  so  that  in 

reality  there  are  only  four  "degrees  of  freedom."  The  four  (inde¬ 
pendent)  coefficients  may  be  used  as  shape  descriptors,  or  all 
seven  may  be  used.  If  d^O  or  d2i0,  however,  the  shape  of  the 

convergence  zone  is  not  suitable  for  description  by  the  form. 
hence,  the  coefficients  calculated  should  not  be  used. 

Therefore,  one  obtains  either  a  set  of  shape  coefficients,  or 
an  indicator  that  the  shape  is  significantly  different  from  the 
anticipated  form. 


r  >  R. 
D 

r  <  R. 
1 


This  completes  the  quantitative  measures  of  a  single  trans¬ 
mission  loss  curve.  The  next  two  chapters  will  deal  with  compari¬ 
sons  of  two  such  curves. 


IX.  Comparison  of  Measures  of  Transmission  Loss  Curves 

If  one  is  given  two  (or  more)  transmission  loss  curves,  a 
natural  way  to  compare  them  would  be  to  compare  the  various  measures 
of  those  curves  which  were  developed  above. 

The  measures  defined  previously  may  be  listed  for  convenience 


w 


l,k 


k,°° 


To  compare  the  monotonic  components  and  &2,  one  might  wish  to 

just  compare  their  coefficients,  but  such  a  procedure  can  be  mis¬ 
leading,  for  the  two  curves  may  be  quite  close  over  the  common 
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interval  [6,  r^2] ,  but  have  rather  different  coefficients.  A 
better  comparison  is  obtained  by  computing 


fr  12  2 

4  (A1(r)-£2(r))^f  (r)dr 


P2 


f (r)dr 


This  integral  is  most  efficiently  computed  by  numeric  quadrature 
rather  than  by  analytic  evaluation. 


This  seems  to  be  a  rather  simple  formulation,  and  taking  the  d 
criminant,  gives 


T6  =  min{  |  Rj_]_“Rj  2  I  :  i,j  =  l,2} 


and  let  iQ,  j  be  the  indices  which  give  this  minimum.  It  will 
also  be  useful  to  form  the  difference  of  the  zone  spacings 


T 


e 


X 


2  * 


These  two  measures  indicate  (by  their  proximity  to  zero)  how  the 
convergence  zone  spacing  of  the  two  curves  compare. 

In  regard  to  the  limiting  values  of  the  ftk(X)'s,  one  might 
be  tempted  to  follow  the  previous  procedures  and  just  take  the 
difference 


But,  since  these  represent  the  level  of  the  first  peak  above  the 
monotonic  component,  the  value  of  this  measure  is  not  immediately 
obvious.  First,  it  should  be  remarked  that  if 

io  *  ^o 

then  no  tests  should  be  performed  on  ^  since  they  correspond  to 
different  peaks  of  the  two  curves.  If  the  indices  are  equal,  then 
form  the  ratio 


Tc  - 


This  is  a  measure  of  how  well  the  first  peaks  correspond  in  light 
of  the  overall  deviations.  If  is  close  to  1,  agreement  is  good. 

Now,  we  consider  the  shapes  of  the  zones.  One  could  compare 
the  curves  zone  by  zone,  but  this  is  probably  more  complicated  than 
the  situation  requires.  It  is  better  to  confine  attention  to  one 
or  two  convergence  zones:  the  first  peak,  the  second  peak,  or  both 


The  choice  of  whether  to  use  the  first  or  second  peak  is  difficult 
since  each  has  certain  advantages,  hence  the  use  of  both.  As  with 
the  previous  test,  this  one  should  only  be  done  if  iQ  =  j  in  the 

calculation.  One  then  uses  a  measure  of  the  form  of  an  integral 
of  weighted  squared  difference  in  the  shapes,  adjusted  for  possible 
differing  levels,  viz. 


where 


o>o  (x)  ('F1  (x+Rj^)  -4*2  (x+Rj2)  +<J>)  2dx 

f  t  A 

J  too(x)dx 


s  =  max  (6-R.,  ) 
k  3K 

t  =  min  (rk-Rjk) 


But,  this  is  rather  unwieldy,  and  if  u>  has  a  reasonably  small  support 
(e.g.,  e<15  nm. ,  cf.  App.  D  and  r12  is  not  very  small  (e.g.,  r12> 

60  nm) ,  then  for  most  transmission  loss  curves,  the  above  functional 
is  equal  to 


T  (j)  =  min 


f*uQ  (X)  (x+Rj  x )  ”4' 2  (x+Rj2)  +♦)  2dx 


/  wQ(x)dx 


For  those  transmission  loss  curves  satisfying  the  above  .assumptions, 
but  where  the  above  functional  is  not  equal  to  the  first  form,  the 
forms  are  usually  so  close  that  the  difference  is  negligible.  The 
evaluation  of  this  measure,  in  closed  form,  is  given  in  Appendix  E. 
The  smaller  this  measure  is,  the  closer  the  shapes  of  the  convergence 
zones  are. 

This  completes  the  measures  based  on  the  measures  of  the  indi¬ 
vidual  transmission  loss  curves. 


X.  Measures  Based  on  Two  Curves 


In  addition  to  the  indicators  based  on  comparison  of  measures 

of  the  individual  curves,  there  are  useful  measures  which  can  only 

be  formed  from  considering  several  transmission  loss  curves  at  once. 

3  6 

For  the  case  of  two  loss  curves,  some  of  these  will  be  given. 

The  obvious  measure,  which  has  already  been  introduced,  is  the 
Lp  norm.  Thus,  we  form  a  weighted  norm 


VP) 


1/p 


and  evaluate  this  for  p=l  and  p=l/2.  The  first  (p=l)  gives  the 
(weighted)  area  between  the  curves,  normalized  by  length  of  the 
common  domain  of  the  curves.  The  use  of  a  value  of  p  less  than  1 
results  in  a  measure  which  is  less  sensitive  to  great  divergences 
(such  as  caustics)  and  quite  sensitive  to  the  small  devaitions 
between  the  curves.  But,  the  integral  (for  p=l/2)  cannot  be 
easily  given  any  physical  meaning. 

To  investigate  shift  in  convergence  zones,  one  might  be  tempted 
to  consider  choosing  x ,  where 

0  *  t  <  r12-  6 

in  order  to  minimize 

(P-l  (r)-P2  (r+x) )  2  (f  (r)  f  (r+x)  )  1//2dr 

[°  (f  (r)f  (r+x))'1/2dr 

'6 


or  to  maximize 


where 


P 


12 


T 


But,  unfortunately,  the  first  form  is  not  nearly  as  good  a  measure 
of  shift  as  is,  and  the  second  form  is  even  worse.  The  second 

form,  which  is  based  on  a  correlation  concept,  is  highly  unsatis¬ 
factory  as,a  shift  measure.  But,  the  perturbation  components  P^fr) 
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of  the  transmission  loss  curves  are  not  periodic  functions.  There 
are  a  number  of  phenomena  that  are  not  present  with  periodic  functions 
which  can  affect  this  form,  in  some  cases  rather  drastically.  Neither 
of  the  two  estimates  for  t  obtained  from  optimizing  these  two  forms 
need  be  reasonable  estimators  of  the  shift.  Thus,  neither  form  will 
be  used;  we  have  discussed  them  merely  due  to  their  rather  exten¬ 
sive  use  in  other  areas  (e.g.,  signal  processing). 


If  the  measure  T  is  not  small  (e.g.,  T  >8  nm)  then  one  may 

fc-  t 

wish  to  examine  the  shape  of  the  convergence  zone — shadow  zone 
cycles  of  the  two  curves,  scaled  so  the  convergence  zones  coalesce. 
Towards  this  end,  let  Vk  be  the  linear  transformation  taking  the 

first  peak  Rlk  into  0  and  the  second  peak  R2k  into  1,  viz 


VR 


lk; 


0 


VR 


2kj 


1 


or 


Vk(r)  =  <r-Rlk)/\ 

These  linear  transformations  take  the  range  onto  a  new  interval, 
on  which  the  peaks  of  the  Pk 1 s  coalesce.  The  new  curves  are  then 
given  by 


and  their  common  domain  [t^,  t2 ]  is  given  by  the  intersection 
[t1,  t2]  =  [V1(6)/  V]_(r1)  ]  H  [V2  (6)  ,  V  2  (r  2 )  ] 

Note  that  since 


t-  >  1 
z 

this  interval  is  not  empty.  Define  the  scaled  weighting  function 

?k(v)  =^7^ 

then  a  useful  measure  of  the  similarity  of  shape  is 

C2  (G1(v)-G2  (v)-c|))2C1(v)  C2  (v)dv 
=  min  - ^ - 

*  /t2  ?1(v)?2(v)dv 

ti 

which  may  be  reduced  to  a  numeric  quadrature  by  proceeding  as  was 
done  in  Appendix  E,  i.e.,  let 

Eq  =  f  2  (G1(v)-G2(v))qc1(vH2(v)dv 

fcl 

for  q=0,  1,  2.  Then 

Tx  =  (E2fi0  -  fi12}/fi02 

The  smaller  is,  the  closer  the  geometric  shape  of  the  convergence 
zones  of  the  two  curves  will  be.  Of  course,  if  any  of  the  requisite 
parameters  (e.g.  ,  ^ ,  X^)  could  not  be  calculated  with  confidence, 

then  this  measure  cannot  be  evaluated. 


Now  one  might  wish  to  compute  characteristics  of  a  pair  of 
transmission  loss  curves  by  first  estimating  £  and  G12  of 

the  curve,  and  then  using  Pkk(r)  rather  than  P^(r)  in  the  calcula¬ 
tion  of  the  various  quantities  in  Chapter  VII  and  VIII.  Such  a 
technique  has  value  when  one  of  the  transmission  loss  curves  is 
defined  over  a  rather  short  interval  (in  comparison  to  the  other) 
and  both  curves  are  a 'priori  believed  to  be  similar.  But,  this  is 
a  refinement  whose  complications  are  not  worth  considering  at  this 
point . 

However,  we  do  mention  one  measure  based  on  2 <  viz 

0k(V 

Vk)  =  7i  "(TT 

®12  '  i2 • 

which  is  a  measure  of  how  well  the  individual  monotonic  components 
compare  with  the  component  ^12*  T^e  sum 


Tk(3)  =  Tk(1)  +  Tk(2) 


is  often  more  handy  to  use  in  practice.  Note  that 


0  *  TK(k)il 


k  =  1,  2,  3  . 

The  closer  (3)  is  to  1,  the  better  the  Jt^2  component  describes 

both  curves,  and  the  closer  their  individual  monotonic  components 
are.  In  a  sense,  T<(3)  measures  much  the  same  characteristic  as 
T  ,  but  there  are  unusual  situations  in  which  one  of  these  measures 

is  numerically  unstable  (i.e.,  a  different  tabular  representation 
of  a  wierd  transmission  loss  curve  could  significantly  alter  the 
measure)  while  the  other  one  will  remain  stable.  Thus,  it  is  pre¬ 
ferable  to  compute  both  Tv  and  the  T  (k)'s. 

T  1C 
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XI.  Implementation  Considerations 

The  next  step  in  a  process  such  as  this  is  to  implement  enough 

of  these  tests  (i.e.,  some  or  all  of  them)  in  order  to  allow  them  to 

be  exercised  against  realistic  transmission  loss  data  of  a  variety 

of  types.  From  such  an  exercise,  meaningful  values  for  the  numerous 

,  2 

parameters  used  in  the  measures  (e.g.,  k  ,  f  (r) ,  Y2,e  ,  d)  should 

be  obtained.  Also,  the  analyst  will  acquire  familiarity  with  the 
different  measures  and  will  see  how  different  types  of  changes  in 
the  loss  curves  will  affect  the  various  measures.  Only  after  such 
a  period  of  analysis  can  one  turn  to  "production  use"  of  such 
measures.  For  a  simplistic  comparison  program,  one  may  quickly 
construct  a  program  functionally  similar  to  Fig.  16,  but  for 
flexibility  and  ease  of  later  conversion  to  an  operational  model,  one 
would  probably  desire  something  along  the  functional  concept  of  Fig. 
17.  In  each  of  these,  the  two  functional  boxes  which  now  concern 
us  are  "Computation  of  Individual  Measures"  and  "Comparison  of 
Curves  (T  measures)":  the  remainder  do  not  lie  in  the  scope  of  this 
paper. 

Under  "Computation  of  Individual  Measures"  one  would  wish  to  com 
pute  *-k/6k»  Xk#  Pk/  #k»  Rlfk'  Ak'  nkfO0  and  ^k-  Now,  the  first  two 

may  be  easily  computed  using  the  determinant  minors  of  Appendix  A, 
where  trapozoidal  integration  is  used  for  the  quadrature.  The  third 
may  be  found  by  testing  on  a  mesh.  Then,  the  decision  as  to  the 
existence  of  convergence  zones  is  made  (cf.  Fig.  22).  if  there  are 
no  convergence  zones,  then  the  remaining  quantities  have  no  meaning, 
and  the  computation  is  complete.  If  there  are  convergence  zones, 
the  process  continues.  The  mesh  used  for  Tk  will  also  be  used  for 

Pk,  and  in  most  cases  can  be  used  for  as  well.  Both  P^  and  4>k 

are  computed  in  tabular  form.  The  next  two,  ^  and  X^  are  local 

maxima  of  functions,  thus  a  common  algorithm  can  be  used  for  both 
of  them.  Of  the  various  numeric  algorithms  available  for  a  local 
maximum  or  minimum,  a  simple  step  search  should  form  a  more  than 


IftSaaflSfr 


adequate  base  for  our  algorithm.  (Due  to  problems  of  multiple 
peaks,  etc.,  several  checks  must  be  added  to  any  "canned"  algorithm 
before  it  is  used  on  these  functions.)  The  measure  ft.  may  be 

found  by  numeric  quadrature,  and  4^  may  be  obtained  from  the  linear 

system  given  in  Appendix  A. 

Under  "Comparison  of  Curves"  the  measures  Tq  through  would 

be  computed.  There  are  ten  types  of  measures,  but  since  some 
measures  are  evaluated  with  different  arguments  (T  ,  Tg,  T^)  there 

are  thirteen  scalar  measures,  and  if  one  includes  the  linearly 
dependent  T^{3),  there  are  fourteen  scalar  measures.  The  first 

seven  (Tq  through  T  )  are  computed  from  measures  of  the  individual 

transmission  loss  curves  (Ch.  IX).  Of  these,  requires  a  quadrature; 

Tg,  T^,  T5,  Te,  and  are  computationally  trivial;  and,  T^  may  be 

evaluated  in  closed  form  using  the  equations  of  App.  E.  The  other 
three,  detailed  in  Chapter  X,  require  numeric  quadrature.  For  TQ  and 

Tl  are  not  integrated  in  closed  form,  while  TK  depends  on  0i2^i2^' 

which  in  turn  is  obtained  from  a  system  of  linear  equations  (App.  A) 
whose  coefficients  were  obtained  by  numeric  quadrature. 

Thus,  it  is  seen  that  implementation  of  these  measures  is  rather 
straightforward.  And,  if  good  use  is  made  of  subroutines  to  handle 
functions  such  as  quadrature,  solution  of  linear  systems,  determinants, 
and  maxima/minima  of  functions,  implementation  should  also  be 
efficient. 

XII.  Conclusions  and  Recommendations 

The  measures  introduced  in  this  report  should  give  a  good 
quantitative  description  of  the  characteristics  discussed  in  Chapter 
III.  However,  there  are  a  large  number  of  parameters  (e.g.,  constants, 
weighting  functions)  which  must  be  determined  by  careful  analysis. 

The  heuristic  values  given  for  some  of  these  are,  as  stated,  merely 
estimates.  Such  an  analysis  would  require  these  measures  to  be 


computed  for  a  wide  selection  of  transmission  loss  curves,  which 
implies  a  computer  program  of  some  kind  to  perform  the  myriad 
calculations.  It  is  strongly  suggested  that  such  a  computer-assisted 
analysis  be  performed  as  the  next  logical  step  towards  an  automated 
procedure  for  determining  similarities  and  differences  of  trans¬ 
mission  loss  curves. 


APPENDIX  A 

Estimation  of  JL 
k 

The  estimation  of  the  minimizing  functions  £^,  ,  and  £12 

follows  the  standard  methods  of  weighted  least  squares  as  developed 
in  Ref.  a  and  extended  in  Ref.  b.  However,  the  functionals  0^,  ©2 , 

and  012  are  formed  by  integrals  of  rational  functions  of  Tk(r),  and 

these  functions  are  only  specified  at  points.  Therefore,  before  pro¬ 
ceeding  with  the  minimization,  it  is  necessary  to  consider  the  effect 
of  this  tabular  data. 

Obviously,  there  is  a  large  class  of  functions  which  will  pass 
through  the  tabulated  points;  thus,  there  can  be  substantial  varia¬ 
tion  in  the  resulting  integrals.  For  most  purposes,  however,  it  may 
be  assumed  that  the  transmission  loss  curve  can  be  modeled  by 

"connecting  the  dots"  of  the  tabulated  data;  i.e.,  linear  inter- 
3  8 

polation.  In  such  a  case,  the  integral  of  T^(r)  is  given  exactly 
by  the  trapozoidal  rule.  That  is,  if  (x . ,  T,  (x . ) )  are  the  range- 

X  K  1 

transmission  loss  pairs,  where 


x  =  6 
o 

xi+i=  xi +  4x 


then 


x  =  r, 
m  k 

m-1 

frk 

1 1  ,  ,  p  ,  1 

/  T.(r)dr  = 

'6  x 

7  Tk(xo>  +  £VX i>  +  7  Vxm> 

1 

Ax 


Now  integrals  of  the  form 


are  not  expressible  exactly  by  the  trapozoidal  rule,  for  f (r)  is 
not  piecewise  linear.  The  error  involved,  however,  in  applying  the 
trapozoidal  rule  is  small,  and  it  is  better  to  use  this  rule  in  such 
instances  than  to  use  either  of  Simpson's  rules. 

With  the  matter  of  quadrature  methodology  addressed,  attention 
turns  to  the  main  problem.  First,  define  four  functions  for  the  loss 
curve  Tk (r ) : 

g*(r)  =  1 
g^(r)  =  log  r 
gk(r)  =  Tk(r) 
g^(r)  =  ar 


and  define  the  determinant  A{k)  of  integrals  as 


Mk)  =  I  ai  j  I  4X4 


where 


k 

aij 


=/; 


k  gk  (r )  g  j  (r )  f  (r )  dr 


Or,  expanding  in  detail, 


lll 


=  [  k  f (r) dr 
A 


a,  _  =  a 


12 

k 


a.  -  =  a 


21 

k 


/rk 

=  J  log  r  f (r)dr 


=  j  Tk(r)f  (r)dr 


*13 

k 


a,  ,  =  a 
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k 
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rrh 

=  pj  r  f  (r  )dr 


k 


=  /  ~k  log2r  f  (r)dr 


=/  K  Tk(r)  log  r  f  (r)dr 

6 


k  C  k 

l42  =  a  /  r  lo9  r  f  (r)dr 


k  rrk  2 
a33=L  Tk  (r)  f  (r)  dr 


k  f  rk 

=  a^3  =  ay  Tk(r)r  f(r)dr 


k  2  r  rk  2  _  ,  ,  , 

a44  =  a  j  r  f(r)dr 


The  coefficients  b,  and  £,  of  £.,  are  given  by  the  system 


Let  M^(k)  be  the  first  minors  of  A(k)39  and  (k)  the  second 
j  mn 
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minors. 

For  notational  convenience,  let 


Nij(k)  =  M34(k) 


Then 


bk  = 


N13(k)+N14(k) 


N34(k) 


bk  * 


N23(k)+N24(k) 


N34(k) 


and  the  minimum  value  of  ®k  is 


29 


APPENDIX  B 

On  the  Spacing  of  Convergence  Zones 


I.  In  order  to  formulate  meaningful  (yet  hopefully  simplistic) 
descriptors  of  the  convergent  zone  placement,  it  is  prudent  to  examine 
a  reasonably  realistic  case  for  which  the  analytical  solution  is 
quite  tractable.  One  may  begin  by  using  standard  ray  theory,  which 
gives  rise  to  caustics  of  infinite  intensity.  While  such  behavior  is 
disastrous  for  propagation  loss  models,  it  is  useful  in  this  situation 
in  that  the  peak  of  the  convergence  zone,  or  the  caustic,  is  easily 
identified.  By  choosing  a  sound  speed  profile  which  is  C  ,  one 
eliminates  the  spurious  caustics  which  sometimes  arise  in  ray  trace 
models  (Ref.  c) . 


II.  Consider  the  analytic  case  with  a  range-independent  parabolic 
profile 


c  ( z )  =  c(zo)/{1-Y2(z-zo)2}1/2 


Y  >  0 


where  the  sound  axis  is  at  depth  zQ  and  no  depth  is  considered 
where 


yIz-ZqI  -  1 


Place  the  source  on  the  sound  axis,  and  for  purposes  of  long-range 
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calculation,  only  consider  the  refracted  paths.  The  ray  paths 

are  given  by  the  equation. 


z  =  z  +  y  sin0^  sin(Y  x  sec©  } 
o  o 


_  i  2 

C  *  Y  cos6^sin{ y  x  sec8_ }  +  x  tan  cos(yx  sec6„) 


where  0Q  is  the  initial  tangent  angle  of  the  ray  (Ref.  c)  ,  and  £ 

is  the  partial  derivative  of  z  with  respect  to  6^  (Ref.  d). 

Normally,  £  is  used  for  the  intensity  calculation.  The  caustics 
will  touch  the  ray  paths  at  the  zeros  of  £;  thus,  the  caustics  will 
be  given  by  the  locus  of  points  (x,  z)  on  the  ray  paths  for  which  £ 
is  zero.  To  simplify  the  analysis,  the  coordinate  system  will  first 
be  translated  so  that  the  origin  is  at  the  source,  viz, 

y  =  z  -  z 
2  o 

Then,  the  locus  of  (x,  y)  will  be  the  caustic.  The  locus  is  given 
by 

y  =  Y  ^sinB^^  sin{yx  sec0Q} 

—  1  2 
0  =  y  cos0^  sin{yx  sec0  }  +  x  tan  0  cos(yx  sec0  } 
o  o  o  o 


The  parameter  @o  is  unwieldy,  so  transform  the  equations  into  a 


new  parameter 


£  =  Y  x  sec0 


to  obtain 


x  =  |  /V1"!  tan^ 


y  =  5ia£  AW 


For  any  value  of  £,  in  the  following  sequence  of  half-open  intervals 
there  corresponds  a  (real)  point  on  a  caustic. 

£  e  ((m-1/2)  tt,  miT] 
in  —  1,  2,  3,  .  .  . 
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The  set  of  caustics  is  plotted  in  Fig.  18.  This  figure  represents 
the  entire  caustic  structure,  which  assumes  that  the  only  boundaries 
are  the  natural  ones  where  c(z)  ceases  to  be  defined  [cf.  the  equa¬ 
tion  defining  c(z)].  If  attention  is  restricted  to  that  range  about 
the  sound  channel  axis  such  that  c(z)  is  within  10%  of  c(z  ),  then 

the  caustic  structure  appears  as  in  Fig.  19,  which  has  been  scaled 
for  visual  purposes. 

Now,  consider  the  spacing  between  successive  caustics  along  a 
horizontal 'cut.  Along  the  sound  channel  axis  the  caustics  have 
cusps  (cf.  Fig.  18)  which  are  separated  by  ir/y.  Along  a  horizontal 
cut  which  is  not  on  the  axis  it  is  seen  from  Fig.  18  that  the 
spacing  is  not  constant,  but  in  the  region  of  interest  to  us  (Fig.  19) 
it  is  nearly  constant.  To  quantify  this,  the  parameteric  equations 
for  the  caustics  are  expanded  about  each  cusp. 


Writing 


£  =  imr  -  n‘ 


one  may  develop  the  power  series 


imr  3  2,3  ' 

X  y  2y^  8miTY  ^ 


1  3 

y  =  — =  n  - 


6y  (imr) 


8m27T2-3  6  ,  8 . 

- =-a-  n  +  O (r»  ) 

4  8m  it  y 


9  11 

n*  +  0  (nxx) 


In  the  vicinity  of  the  cusp,  the  branch  of  the  curve  may  be  written 
(Bk.  II,  Ch.  II  of  Ref.  e)  as 

x  =  f  -3  ^y2/3  +  3  ^Ty4/3+0(y2) 


Let  T  denote  the  constant 


and  consider  the  spacing  between  caustics  of  the  transmission  loss 
curve  at  constant  depth 


•.*  ■_«  V-  \ 


Figure  18.  Caustics  for  Parabolic  Profile,  Natural  Boundaries 


for  u  small  in  comparison  with  both  the  total  depth  and  1/y.  Then, 
the  spacing  between  the  m-th  and  (m+l)st  caustic  at  that  depth  is 
given  by 


+  0 (u2 ) 


which  may  be  approximated  by 

V  .  r3  -  \  r  m'2'3  u2/3  -  ir'1  m'4'3  u  4/3 

But,  since 

M  <<  i/y 

it  is  clear  that  A  x  is  almost  a  constant. 

m 

III.  Now  consider  a  case  where  the  source  does  not  lie  on  the 
sound  channel  axis.  Let  the  sound  speed  profile  be  the  same  as  above, 
but  place  the  source  at  range  0  and  depth  z^  where 

Z1  *  zo 

Let  be  the  tangent  angle  of  the  ray  at  the  source,  and  let  ©  be 

o 

the  tangent  angle  of  the  ray  path  at  the  sound  channel  axis,  so  they 
are  the  ray  angles  at  depths  z1  and  zQ,  respectively.  Then 


c(z,)  sec6..  =  c(z  )  sec0 


The  ray  paths  are  given  by 


z  =  zQ  +  Y  1  jEsinfy  x  sec9o}  -  yUq-z.^)  cos{y  x  secBo}( 

where 

2.2  2  2 
5  =  sin  6q  -  y  (zQ-z1) 

so  that 


sin6  cos6 
o  o 


and 


- „  sinQ  cosB 

Yyg-  =  - = -  sinfy  x  sec6o)  + 

o  ~ 

H  y  x  tan6^  sec©^  cos(y  x  secB  }  + 
o  o  o 

Y(z  -z,)yx  tan0„  sec0^  sin{y  x  sec6  } 
O  J.  o  o  o 

Now  d0o/d0^  is  non-vanishing  for 

lzl"zol  <  1/Y 

,  42 

so  instead  of  setting  c  to  zero,  one  may  set 


Making  the  transformation  to  the  more  tractable  parameterization 

£  =  Y  x  sec6 


which  was  used  in  the  previous  case,  yields 


:2  =  1  -  yV  -,Vr2 


Since 
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sin6  ^  0 
o  ' 


it  may  be  divided  out  of  the  equation  for  the  partial  derivative, 
which  then  becomes  (upon  multiplication  by  xH) : 


0  =  x  cose  sin{yx  sec  }  + 
o  o 


2  2  2 

(z--z,)y  x  Esec  6  sin{y x  sec6  }  + 
U  1  o  o 


2  2  2 

yx  E  sec  6  cos  {yxsec6  } 
o  o 


or,  in  the  new  parameters, 


sin?  +  u?2Esin?  +  y_1?2E2cos?  =  0 


For  notational  simplicity,  a  number  of  Bj's  will  be  defined,  starting 
with 


Bl  =  -yV2 


2  2 

B2  =  l-y2u2 


Rationalizing  the  equation  for  x  gives 


(y4sin2?)x4  +  2y2?3sin?cos?  (B^+BjX2)  +  ?6cos2?  (B12x4+2B1B2x2+B22) 


=  y2u2?6sin2?  ( B-,  x2+B~) 


which  may  be  written  in  powers  of  x  as 


8,  +  B.x2  +  Bcx4  =  0 
3  4  5 

where 

83  =  56622cos2S  -  Y2u2£6B2sin2£ 

B4  =  2y2£382Sin£cos£  +  2B^82^cos2£  -  y2u2 8-^sin2 £ 

85  =  y4sin2C  +  281y2£3sinCcos£  +  812C6cos2C 
or,  in  their  factored  form 

63  =  Is  (l-y2u2) (cos2?  -  y2u2) 

64  =  y2S3{2cosUsinS  -  £cos£)  (1-  y2u2)  +  y2u2£sin20 

85  =  y4(sin£  -  £cos£)2 
Now,  make  the  change  of  variable 

*  =  xVr3 

so  the  polynomial  transforms  into  the  quadratic 

Bg  +  B7^  +  6g4>2  =  0 

where 

Bg  =  (l-y2u2) (cos 2 £  -  y2u2) 

8?  =  2cos^(sinC  -  £cos£) (l-y2u2)  +  y2u2£sin2£ 


2 


(sin£  -  £cos£) 


This  seems  to  be  a  rather  simple  formulation,  and  taking  the  dis¬ 
criminant,  gives 


e72  -  466B8  =  {2yusin£}2B92 


where 


B92  =  (l-y2u2  +  i-y2u2£2}sin2C  -  £  ( l-y2u2 )  cos£sin£ 
The  roots  of  the  quadratic  are 


* 


-B7±2YuB9sin£; 
2  8  8 


Thus,  to  compute  the  caustics,  one  uses  the  parameter  £,  computes 

ip,  then  E,  then  x  and  y,  remembering  that  most  of  the  functions  are 
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multivalued,  i.e., 

E  =  V  1-Y2u2-^ 
x  =  y"1^/!? 

y  =  y_1  { Esin£~YUcos£ } 


Now,  for  £  a  positive  integral  multiple  of  it,  there  is  a  (real  ) 
cusp  on  the  caustic,  which  may  be  obtained  by  power  series  develop¬ 
ments  Let 


then 


£  =  IMT 


meZ 


+ 


o(n 


16 


) 
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V  s 


4  1 


By  =  -2 7Tm82  +  71111  (2-y"u“)  n  -  ■j(2+y-'J")  n”  -  j  nm (2-y-u“ )  n~  + 


1  ,,  2  2.  10  ,  2THT1,-  2  2,  12  14, 

^5(6-7  u  ) n  +  T5“^2-y  u  ^ n  +  0  n  ^ 


„  _2  2  _2  2  4  2TTm  6  .  tt  m  8 

Bg  =  'fTm  -  Trrnn  +  — j-  n  +  — ^ —  n 


2  10  ,  1  ,,  ,  2  2, 
g-  7rmri  +  ^-(5-2tt  m  )n 


+  o(n14) 


0  ,,  2  2,  2,1  2  2  2  ,2  4  .  _  ,  6, 

Bg  =  TTm(l-Y  u  )n  +  j  y  inii  M  +  0(n  > 


tJj  =  —  + 
irm 


!i  2yuV^2 

™  +  /„^3/2 


n3  -  n4  +  0 (n5) 


_  2  2  3TrmB-j  0  2  (nm)  3/2ix/b7 

x  - - —  e2 - 2~  n  + - Y - n 

Y  Y 


+  {3y“282  -  |  m27T2u2}n4  +  0(n5) 


so  the  development  of  the  curve  about  the  cusp  is 


x  -£wjin3 

j  =  0 


y 

3=0 


where 


w0i  ■  kitty 


t  i »m+l 

w02  *  (-1)  u 


wu  - 0 


39 


0 


w 


12 


w 


21 


-  1  y-1 
2  Y 


V®7 


w22  “  0 


w 


31 


=  u  ^imt 


w 


32 


w 


-  2  2  2  2 
3  B2  “  2m  tt  y  o 


41 


w 


42 


8ym7 


(-l)m  j 


and  so  forth. 


Thus,  these  cusps  lie  at  uniform  intervals  in  the  range 
coordinate,  but  alternate  in  depths  of  y=±u.  At  one  of  the  depths, 
there  are  cusps  at  a  uniform  spacing,  but  the  caustic  branches 
between  the  cusps  are  not  equally  spaced.  But,  since 

M  <<  i/y 


and 

e2  -  i 

the  spacing  is  close  to  uniform. 

IV.  The  final  case  to  be  considered  is  more  realistic  in  that  it 
has  a  surface  boundary.  Let  the  profile  be  parabolic,  as  before, 
with  the  channel  axis  at  depth  zQ ,  where 

0  <  zq  <  1/y 
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As  before,  the  depth  coordinate  is  changed  to 


Y  =  z  -  z. 


so  that  the  source  is  at  the  origin,  and  let  0q  be  the  tangent 


angle  at  the  source.  Let 


T.  =  y  cos0^  Arcsin(yz  cos0  ) 
1  o  '  o  o 


r2  =  +  y  ttcos0o 


-  ■  ra 


eio< 


where  [  ]  is  the  greatest  integer  function.  And,  define  the  threshold 

=  Arcsin(yz0) 

Then,  the  ray  paths  have  two  forms,  depending  on  whether  they  are 
refracted  or  surface  reflected  and  refracted: 

y  =  y~1sin8^sin{y x  sec6  }  I  6  I  £  F 
o  o  o  o 

y  =  y-1sin6osin{y (x-B10 (x) r2) seceQ}  Fo<|0o|  <  J 

Thus,  the  locus  of  caustics  is  described  by  two  sets  of  equations, 
corresponding  to  the  two  cases  above.  In  the  first  case,  the  caustic 
may  be  described  by  the  same  parametric  equation  used  in  the  first 
case,  except  the  domain  of  the  parameter  £  is  different. 

In  the  previous  case,  the  domain  was 

£e  ((m-l/2)fr,  imt]  meZ+ 

But,  in  this  case,  the  parametric  equation  for  x  yields 


1  i  ’  O';  O  *  * 


cos  0 

o 


or 


£cot£  +  cot  e 


0 


or 


E,cotE,  +  cot28Q  £  0 

Let  be  the  root  of  this  equation  between  (m-l/2)Tr  and  mir . 
in  practice 


zo  <<  1/y 


the  solutions  are  approximated  by 


=  mTr{l-tan2eo  +  tan46Q  +  (|  m2ir2  -  1)  tan60o 


+  (l-^m2Tr2)  tan80  -  (im47i4  -  m2ir2  +  1)  tan106  } 
J  O  3  J  O 


+  0 (tan140  ) 
o 


And,  the  domain  of  the  parameter  is  then 


Se(£m,mTT]  mez+ 

Thus,  along  the  sound  channel  axis  are  the  cusps  shown  in  Fig.  18 

Now,  consider  the  caustic  of  the  second  part  of  the  ray  path 
family.  Since  interest  will  be  restricted  to  the  vicinity  of  the 
sound  channel  axis,  it  is  sufficient  to  just  consider  those  por¬ 
tions  of  the  caustic  which  touch  or  cross  the  axis.  This  will 
require 


•V 
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simultaneously.  But,  from  the  second  form,  it  is  obvious  that  for 


0Q  ^  0  modiT 

these  two  cannot  both  be  satisfied.  Hence,  along  the  sound  channel 
axis,  there  is  a  series  of  equally  spaced  cusps  on  the  caustic. 
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APPENDIX  T 
On  the  Choice  of  w 


The  short-span  smoothing  function  used  in  the  calculation 

of  ^(R)  must  be  defined  before  numeric  compuations  can  be  performed 

Fortunately,  there  is  a  very  wide  latitude  in  the  choice  of  such  a 
function,  and  the  final  results  are  not  very  sensitive  to  the  choice 
of  function.  Theoretically,  one  may  choose  eoR(r)  to  be  almost  any 

C0“  function  whose  support  was  confined  to  a  sufficiently  small 

interval  about  R,  say  [R-e,  R+e]  .  If  uu.(r)  is  C  °°,  then  $>  (R)  is  C°° 

K  O  K 

and  since 

supp  [R-e,  R+e] 

00  00 

$k(R)  is  CQ  .  In  practice,  however,  CQ  functions  generally  do  not 

lend  themselves  to  rapid  computation  and  one  really  does  not  need 

00  2 
$k(R)  to  be  CQ  ,  but  merely  C  . 

Thus,  two  possible  forms  for  wR(r)  are  given,  one  reasonable 
(but  not  too  good)  and  one  excellent. 

Form  Is  Choose  constants  y ^  where  is  a  dimensionless  quantity 

greater  than  1  and  y 2  is  a  positive  distance  (in  miles) ,  such  that 

e  =  VY1"1)/Y1 

is  a  reasonable  distance  for  integration.  Then  define 

Y7  = 


( Y 1  “  1 )  /  Y 


i.3 

Ml 


"R(r)  “ 


Y1X+(1-V1)  if  x  c  y- 


if  X  i  Y' 


Y2 (x-1) ' 


Vr)  =1-7 - ~ 

(x-1 ) 


if  x  <  and  R  £  r 


if  x  <  y3  and  R  >  r 


2  3 

y2  (x-ir 


*R(r)  — 


if  x  >  Y- 


if  x  <  y3  and  R  £  r 


if  x  <  y3  and  R  >  r 


Yj  (x-1) 
0 


if  x  i  Y' 


Then,  the  supports  of  the  functions  are  intervals  of  width  2e, 


supp 

“r 

[R-e, 

R+e) 

supp 

UR 

[R-e, 

R+e] 

supp 

n 

Pi 

m3 

[R-e, 

R+e] 

Also 


3w_ (r) 


3  wR(r) 

3R2 


“R(r) 


except  at  the  three  points 


r  =  R  -  e 
r  =  R 
r  =  R  +  e 


where  the  left-  and  right-hand  first  partials  exist,  but  are  not 
equal,  and  no  second  partials  exist. 

A  reasonable  set  of  values  is 


Yl  =  2 


which  gives 


Then 


Y2  =  8  nm. 


c  =  4  nm. 


x  =  8  |R'r 


u)R(r)  =  max  ( 0 ,  2x~j~) 


d  d$k(R)  fr} 


KM  =  —  =  /  k  Pk(r)uR(r)f  (r)dr 

6 


k'"  dR 


n  d  d  $k  (R) 


VR)  *  — 2 

K  dR 


=  /  ~  Pk(r)SR(r)f  (r)dr 


+  J{Pk(R-4)f(R-4)Ij6frk] (R-4) 


+Pk(R+4)f (R+4)I [6  r  j  (R+4) } 

JC 


-  &  Pk(R)f(R)IIdfrk](R) 
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where  I  is  the  indicator  function,  viz 

0,  if  r  t  A 

IA<r>  - 

1,  if  r  e  A 


2 

Now,  since  P  is  continuous  (though  not  differentiable),  4>k(R)  is  C  , 

but  the  rather  jagged  form  of  Tk(r)  may  give  (R)  a  very  large  Total 

Variation  which  has  some  undesirable,  but  not  fatal,  side  effects 
in  terms  of  computer  implementation. 


Form  II.  Choose  a  constant  y ^  to  be  a  positive  distance,  in  nm. ,  such 
that  e  =  ny^  is  a  small  number  of  miles  (e.g.,  4  nm. ) .  Define 


x  = 


R-r 


/  1  +  cos  x 


WrU) 


=  < 


x  <  TT 


X  >  TT 


“R(r) 


sin  x 


X  <  IT 


X  >  TT 


“R(r)  = 


COS  X  X  <  TT 


X  >  TT 


Then 


3wR(r) 

5r 


'  “R(r> 


at  all  points 


3  ojr  (r ) 


=  u(r) 


except  at  r  =  R  ±  e,  where  left-  and  right-hand  second  partials  exist, 
but  are  not  equal.  In  this  case,  we  have 


*£(R)  =)  Pk(r)SR(r)f  (r)dr 
6 


k (R)  =/rk  pk(r)SR(r)f(r)dr 


This  form  for  to  gives  rise  to  a  much  smoother  $£(R)  (indeed,  it  is  even 
C1) ,  which  will  not  have  the  problems  that  can  arise  with  the  previous 


Appendix  A 

Estimation  of  4>.  . 

k 


3 


The  shape  of  the  transmission  loss  curve  T,  about  R.  was  modelled 

K  3 


using  the  form 


(r)  H 


dl  +  d2X  +  5^+d" 


ld4  +  d5x  + 


r  i  R. 
3 


r  <  Rj 


Since  the  monotonic  component  should  not  affect  the  coefficients  of 
the  shape  function,  ^(r)  will  be  fit  to  the  perturbation  rather  than 

the  loss  curve  itself.  But  Pk(r)  may  be  quite  rough  (cf.  Fig.  8);  thus, 
rather  than  fit  the  smooth  curve  (r)  to  Pk  (r) ,  it  will  be  fit  to 
$k(r),  the  smoothed  perturbation  function. 

The  form  for  was  chosen  so  that  a  variety  of  shapes  common 

to  convergence  zone  peaks  could  be  accommodated.  In  this,  however, 
the  parameter  dQ  should  not  be  estimated,  but  should  be  chosen  before¬ 
hand.  This  ensures  that  when  one  is  comparing  different  transmission 
loss  curves  T^  and  T2»  then  the  examination  of  the  convergence  zone 
shapes,  and  can  be  easily  accomplished.  For  if  dQ  were 

different  for  and  meaningful  comparison  would  be  very  difficult. 
Thus,  dQ  will  be  a  chosen  constant,  say  5  nm.  (this  may  be  adjusted 

when  implemented).  Furthermore,  should  be  C1,  but  not  necessarily 

2  K 
C  .  Thus,  one  has 


d,  +  d3/d0  =  dA  -  dc/d 


6  o 


d2  "  d3/do2  =  d5  ‘  d6/do2 


The  first  equation  is  the  condition  that  ^(r)  he  continuous  at 

R  ■  ,  the  second  equation  is  the  condition  that  H'/  (r)  be  continuous 
J  k 


at  R j .  These  may  be  rewritten  as 


d3  ~  d6  =  do  (d2’d5) 


d3  +  d6  “  do^d4-dl^ 


Now,  let  uL(r)  be  a  weighting  function  (as  developed  in  App.  T) , 


which  need  not  be  identical  with  uj  (r)  .  This  new  weighting  function 


is  to  allow  the  fit  to  be  made  in  the  vicinity  of  the  convergenze 
zone,  and  allows  for  situations  such  as  Fig.  20  where  the  zone  shapes 
are  the  same,  but  their  size  is  different.  For  this  type  of  fit,  a 
weighting  function  of  Form  I  (cf.  App.  D  is  superior  to  one  of 


Form  II. 


We  then  wish  to  choose  the  coefficient  of  to  minimize  the 


form 


n  =/ 

'X 


l6  uRj(r)|$k(r)-'fk(r)]‘dr 


This  may  be  quickly  done  using  Lagrangian  multipliers,  by  forming 


2  n  +  d7  [d3+dg+d0  (d1-d4)  ]  +  dg  [d3-d6-dQ  (d2~d5)] 


where  d7  and  dg  are  the  multipliers.  Then,  define  the  functions 


{r-Rj+do) 


Vr) 


and  the  integrals 


Jpq  =J6  WRj(r)epeqdr 


cu^  (r)e  e  dr 

Rj  p  q 


Define  the  column  vectors 

.  fd.7 


and  the  8x8  symmetric  matrix 


11 

J12 

J13 

21 

J22 

J23 

31 

J32 

J33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

511 

512 

513 

521 

^22 

^23 

531 

532 

533 

"do 

0 

1 

Then,  the  vector  D  which  minimizes  the  form  is  given  by  the  linear 
system 

HD  =  F  . 
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By  inspection,  one  may  reduce  this  to  a  4x4  system,  but  the  previous 
form  has  advantages  (it  is  a  symmetric  bordered  block  diagonal  matrix) 
For  this  curve  fit  to  be  geometrically  meaningful,  one  should  have 

d2  <  0 
d5  >  0 

If  either  of  these  conditions  is  not  met  by  the  solution  vector  D, 
then  the  convergence  zone  has  a  shape  which  is  not  amenable  to  this 
type  of  fit,  and  the  D  vector  should  not  be  used. 


Appendix  E 
Evaluation  of  T 


In  this  appendix,  the  problem  of  computing 


T  ( j) -=  min 


.  00 

L"  w0(x)  (f1  (x+Rjl)-¥2  (x+Ri2)+$)2dx 
00 

I  <jjq  (x) dx 


will  be  considered.  For  notational  convenience,  let 


JO 

Eg  =  /  “0(X)  (H'1(x+R,1)-Y2(x+R,2))qdx 
E"  =  /°SQ  (x)  ('l'1(x+Rjl)-'i'2(x+Rj2))qdx 


E  =  E  +  E 

q  q  q 


q  —  0 ,  1 ,  2 


Using  standard  techniques,  the  minimum  occurs  when 


*  =  -E1/Eq 


and  is 


Vj)  =  (E2E0  -  E^l/Ej,2 


where  the  E's  are,  of  course,  dependent  upon  the  index  j 


A. 


<  ,  « 

’o 


E1  Y2Y3{h4Yl"2Y2 (Yl+1)h5} 


dTTT  {Y1d0+Y1Y2-Y2}{^n(Y2Y3+do)-£n  dQ} 
O  2 


$ 


g-|—  {(-h4+>2h5)  (d0+Y2)+h6}tnVl 

O  2 


j 

*kr« 


-£  ~  2  Y  -j  hih,  2  2  ^Y1^2^3X 

E2  =/  ^Ylhl  +2Y1h1h2x  +  +  Yxh2  X  + 

o  o  o 

^y2^1^3 


Ylh3 

(x+do) 


7  + 


Y  2^l2  2Y  2^1^2X 


x-y2  x-y2  (x+do)  (x-y2) 


^Y2^2^3X 


y2h3 


+  Y 2^2  x— v  ^  (x+d  V  f x-Y  Y"  ^  2 

2  2  2  x  y2  (x+dQ)  lx  Y2)  (x+d  )  yx-Y_) 


}dx 


2  2  2 

=  y^Y2Y2  +  Y2Y2  Y3  ^2h2  +  2Ylhl h3^ 8,11  ^Y2Y3+do^ -S,n 

+  j  Y^Y2^Y33h2^  +  2Y^^2h3 { Y2Y3~d0  to  (Y2Y3+d0)  +  dQ  to  dQ} 


,  £ 

Y  y  y  n  *y 

+  d1  (YoYo+d  )  ~  Y2hl  Zn  Y1  +  2y2  hlh2{Y3-JlnYl} 

O  fc  J  o 

♦  >2!do?  <*■>  ao  '  in^21,3+ao)  •  myl5  +  Y23h22{lY32+Y3-tnYl) 
+  <Vn(Y2Y3+'Jo>  -  do  ta  do  ’  y2  ta  V 

2  O 


Y2*V 


<Vdo' 


do+Y2  d0+Y2 


VY2Y3 


+  to  dQ  -  to  (Y2Y3+d0) -toY1j 


(Yl”1)Y2 


hl2+Y2Y3hlh2  +  3  Y2  Y3  h2  +  2h2h3  +  d„  (YoY^+d,, ) 


‘o' T  2  T  3  o 
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If  wc(r)  has  Form  I  of  Appendix  T (which  is  the  preferred  form 
for  this  function,  cf.  App.F),  then  one  has  (reiterating  App.F) 


Vx>  =  Yi  +  i/{-i+M/y2> 

y3  =  1  -  1/Y, 

£  =  Y2Y3  >  0 


Let  h^  be  the  difference  of  the  d^'s  for  the  j-th  peak  two  transmission 
loss  curves,  i.e.. 


hi  =  di(1)  "  di(2) 


i  =  1,  2,  3,  4,  5,  6 


where  d.  (k)  is  d.of  T,  (r)  at  R . . 

1  1  K  J 


The  evaluation  of  the  integrals  is  straight-forward. 


E0  *jf(Vj%)ax  "  Y1Y2Y3  -  VnYl 


E0  ■  ES 


^1^3  Y2^1  y2^2x 


Y2h3 


(Y1h1+Y2h2x  +  +  x!./  +  (x+dj  U-y.)  }dx 


1  2  2 

=  Y3Y2Y3h-^  +  j  yiy2  y3  ^2  +  Y^h3  {  (Y2Y3+dQ) -&n  do^ 

2  Y  2^3 

-  Y2h1JlnY1  +  Y2  h2{Y3-^nY1>  +  g  +y^  Un  dQ~ In  (Y2Y3+dQ)  -toY1  > 


=  Y2Y3{h3Y3  +  ^  Y2(Y1+l)h2> 


+  d“+77  {Y1d0+YlY2-Y2}Un(Y2Y3+do)-iln  do} 
O  2 


~  -T 


Y  2 


{  (h,  +Y0ho)  (d  +Y,)+h,}JtnY, 


+  Y2  Y3  2  Y2(Y3+2)h2  +  2h1h2  -  dQ  (Y2+d0)  (Y2Y3+dQ) 


+  {2(Y2(Y1-1)  +  Yxd0)  (hfdoh2)  (Y2+dQ)  -  Y2h3) 


(W 


•  Un(Y2V<V  -  £n  do} 


h-  2 

Y2  hl  +  Y2h2  +  y2+do  £nYl 


®2  ~  ^Y1~^^Y2  ^4  ~  Y2Y3^4^5  +  3  Y2  Y3  ^5  +  2^5^6  +  d  (Y2Y3+d  ) 


+  Y22Y3  j  Y2(Y3+2)h52  -  2h4h5  -  do(Y2+d^) (Y2Y3+dQ) 


-  {2(y2(y1-D  +  Yxd0)  (h4+dch5)  (Y2+dQ)  +  Y2h6> 


{£n(Y2Y3  +  dQ)  -  in  dQ} 


(Y2+do} 


l  h6  2 

Y2  |h4  Y2h5  "  Y2+dQ  £nYl 


From  these  terms  T  (j)  may  easily  be  computed. 


Notation 


Since  a  number  of  different  mathematical  symbols  are  used 
in  this  paper,  this  list  of  notation  and  variables  is  given  for 
convenient  reference. 

d  is  defined  to  be  equal  to 

=  is  identically  equal  to 


Latin  Notation 

A(k)  determinent  used  to  obtain  monotonic  component  of  T^, 
defined  in  App.  A. 

V 

element  of  A(k) ,  defined  in  App.  A 
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coefficients  of  defined  in  App.  A 


spaces  of  continuous  curves 
sound  speed  profile,  used  in  App.  B 
column  vector  used  in  App.  A 
coefficients  of  shape  function  ^ 


.integrals  of  polynomials  in  the  defined  in  App.  E 

K  r 

integrals  of  polynomials  in  the  G.  functions,  defined 
in  Ch.X.  -1 


functions  used  as  integrand  factors  in  App.  A 


column  vector  of  integrals,  used  in  determination  of 
V,  ,  defined  in  App.  A 

weighting  function  for  integrals,  defined  in  Chap.  IV 
scaled  form  of  the  perturbation  function,  defined  in  Ch.X 


functions  used  as  factors  in  integrals,  defined  in  App.  A 

symmetric  matrix  of  linear  system,  defined  in  App.  A 

differences  in  the  coefficients  of  the  two  ¥  's,  defined 
in  App .  E . 

indicator  function  for  set  A  with  argument  r,  used  in 
App.  r 

index  optimizing  T^,  defined  in  Ch.  IX 


integrals  used  in  the  matrix  H  and  vector  F,  defined 
in  App.  A 


index  optimizing  T^,  defined  in  Ch.  IX 

index  of  transmission  loss  function 
Lp  norm  of  function,  defined  in  Ch.  II 

montonic  component  of  loss  function,  cf.  App.  A 
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(k) 

Mil  (k) 


N. . (k) 
13 


Qx(r) 


R1 ,  R1 ,  k 


Tk  (r) 


first  order  minor  of  A(k) ,  used  in  App.  A. 

second  order  minor  of  A(k) ,  used  in  App.  A. 

index  of  caustic,  used  in  App.  B. 
abbreviated  minor  of  A(k) ,  defined  in  App.  A. 

perturbation  component  of  Tk,  defined  in  Ch.  IV. 

perturbation  component  of  T,  with  respect  to  i,--. 

K  12 

norm  value  used  for  TQ 
parameter  used  in  defining 

K 

function  used  to  determine  spacing  of  sucessive  convergence 
zones,  defined  in  Ch.  VII. 

parameter  used  in  locating  first  convergence  zone 

distance  from  source  to  first  convergence  zone  (actually, 
to  the  first  peak  of  P,  which  may  or  may  not  be  the  first 
convergence  zone) 

distance  from  source  to  jth  peak  of  Pk 

maximum  range  of  tabulated  data  for  kth  transmission 
loss  curve 

maximum  common  range  of  transmission  loss  curves 

lower  limit  of  integration,  used  in  Ch.  IX 
k  transmission  loss  function 


t  upper  limit  of  integration,  used  in  Ch.  IX 

I^i  tj]  intersections  of  domains  in  uniform  range  space,  defined 
in  Ch.X 

Vk  linear  transformation  mapping  first  two  peaks  of  p,  into 

0  to  1  respectively,  defined  in  Ch.  X  k 

v  argument  in  uniform  range  space,  used  in  Ch.  X 

wij  Taylor  series  expansion  coefficients,  defined  in  App.  B 


■IV.' vV.. •.-S' 


ro  i 


x 


general  variable  used  for  range,  particularly  for  caustic 
calculation,  cf.  App.  B 

y  depth,  measured  from  sound  channel  axis,  used  in  App.  B 

Z+  the  positive  integers  ^1,2,3,..^ 

2,2  ,z  depths  for  use  with  sound  speed  profiles,  cf.  App.  B 


Greek  Notation 


A  . 
D 


attenuation  coefficient 

variables  of  convience,  used  in  second  case  in  App.  B 
constant  used  in  App.  B 

coefficients  used  in  representations  of  to  and  to  functions, 
defined  in  App.  T 

parameter  of  parabolic  sound  speed  profile,  used  in 
App.  B 

foward  difference  operator 

minimum  range  of  tabulated  transmission  loss  curves, 
defined  in  Cii.  IV 


e  parameter  related  to  support  of  to,  defined  in  App.  r 

perturbation  component  for  convergence  zone  spacing 

5  derivative  whose  zeros  identify  caustics,  defined  in 

App .  B . 

(v)  scaled  weighting  function  for  uniform  range  space,  defined 

K  in  Ch.  X 

H  inf initesmal  used  in  development  of  power  series  for 

caustics,  defined  in  App.  B 

06  minimization  functionals  for  monotonic  components, 

k  1  defined  in  Ch.  V 


0 

K,K 


angle  of  raypath  at  depths  z  and  z1  for  caustic 
calculations,  used  in  App  B  °  J" 

threshold  for  xk>  used  in  Ch.  VI 


thresholds  for  tests  dealing  with  existance  of  convergence 
zones,  used  in  Ch.  VI 


••  ■  A  . •  •- 


T  ,etc 
a 


'average'  distance  between  convergence  zones 
estimate  of  A  for  transmission  loss  curve 


geometric  mean,  used  in  Ch.  IV 

form  used  in  raypath  and  caustic  calculations  for 
second  case  in  App.  B 

used  in  parameterization  of  caustic  curves,  defined 
in  App .  B . 

roots  of  a  transendental  equation  used  for  caustics, 
defined  in  App.  B 

optimization  functional  for  4^,  defined  in  App.  A 
constant 

upper  limit  of  integration  for  a  candidate  measure, 
used  in  Ch.  X 

threshold  for  cr^,  used  in  Ch.  VI 

optimizing  value  of  parameter  a  for  x^/  defined  in  Ch.  V 

shift  parameter  in  candidate  measure,  used  in  Ch.  X 

measures  relating  to  the  similarities  of  the  two  trans¬ 
mission  loss  curves 

parameter  for  depth,  used  in  App.  B 

smoothed  perturbation  function 

minimization  variable  used  in  functionals 

function  used  for  describing  shape  of  peak,  defined 
in  Ch.  VIII 

transformation  variable  used  in  App.  B 
measure  of  variation,  defined  in  Ch.  V 

function  used  for  determining  convergence  zone  spacing, 
defined  in  Ch.  VII 
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limiting  value  of  Q 


weighting 


functions 
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A  metric  is  a  generalization  of  the  concept  of  distance. 

It  is  a  function  of  two  elements  of  the  space  (in  this 
case  transmission  loss  curves)  into  the  nonnegative  real 
numbers.  A  number  of  metrics  exist  for  function  spaces, 
but  it  is  extremely  difficult  to  give  a  meaning  to  the 
metric,  as  it  applies  to  transmission  loss  curves. 

The  problem  of  comparing  a  number  of  transmission  loss  curves 
is  not  significantly  different  from  that  of  comparing  two 
such  curves.  One  may  compare  the  curves  pairwise  and  then 
form  composite  measures  from  all  possible  pairs.  Or  one  may 
select  one  curve  as  the  'base  curve'  and  compare  all  other 
curves  to  it.  The  details  of  the  multiple  curve  comparisons 
will  not  be  addressed  in  this  paper. 

The  need  for  an  automated  procedure  precludes  having  an 
acoustician  within  the  processing  loop  making  somewhat  subjective 
judgements.  But  this  in  no  way  implies  that  the  judgment 
of  the  analyst  is  not  needed.  Critical  evaluation  of  the 
results  by  the  analyst  is  a  requirement,  and  these  measures 
serve  as  a  tool  to  indicate  what  the  types  of  differences 
are,  and  how  great  they  are..  But  these  measures  will  not,  by 
themselves,  ever  explain  why  a  particular  characteristic  is 
present  or  absent.  The  requirement  for  reasonably  efficient 
methods  of  computation  rule  out  some  of  the  more  esoteric 
methods  of  pattern  recognition. 

The  problem  of  comparing  two  transmission  loss  fields  is  not 
addressed  to  this  paper.  If  the  two  fields  were  grids  whose 
set  of  depths  matched,  then  one  could  compare  the  fields  by 
comparing  the  two  transmission  loss  curves  at  each  common 
depth. 


5. 


To  circumvent  this  problem,  it  will  be  assumed  throughout  this 
paper  that  the  transmission  loss  curve  is  a  linear  spline  whose 
grid  points  are  the  tabulated  values  of  the  transmission  loss. 
In  other  words,  the  transmission  loss  curve  is  obtained  from 
the  transmission  loss  points  by  'connecting  the  dots'.  More 
sophisticated  interpolatory  techniques  could  be  utilized  in 
place  of  straight  line  segments,  but  the  advantage  appears  to 
be  minimal  and  the  resulting  complexity  (in  both  interpolation 
and  quadrature)  is  substantial. 


6.  i.e.,  for  any  r  in  the  interval  [a,b]  both  T^r)  and  T2(r) 

will  be  defined. 


7. 


L  forms  a  norm  only  for  p>l,  but  the  form  (and  the  function 
P  ~  '* 

space J  is  defined  for  any  p>0.  For  p  between  0  and  1  the 

space  is  not  a  Banach  space.  Note  that  and  correspond 

to  standard  deviation  and  maximum  deviation,  which  form  the 

basic  measures  of  perturbation. 


8.  Thus,  any  concept  of  'closeness'  for  transmission  loss  curves 
must  take  into  account  the  phenomena  of  shifting  caustics 
such  as  in  Figure  1. 


9.  A  point  criteria  is  any  sort  of  measure  or  functional  which 

depends  on  examining  the  functions  T^(r)  at  fixed  range  points, 
rather  than  viewing  them  as  elements  of  a  function  space. 

Since  the  L  norms  measure  differences  in  the  T,  (r)  at  various 

P  K 

points  (albeit  an  infinite  number  of  points) ,  they  are 
essentially  point  criteria. 

10.  Obviously  the  range  to  which  a  transmission  loss  function  is 
known  is  a  major  concern,  and  if  it  is  too  short  the  true 
nature  of  the  curve  may  not  be  apparant.  But  a  scheme  of 
analysis  must  not  be  so  unstable  that  the  results  are  highly 
dependent  upon  whether  the  transmission  loss  data  cuts  off 
at  a  convergence  zone  or  a  shadow  zone. 
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11.  If  one  is  dealing  with  shallow  water  propagation,  RAP  propa¬ 
gation,  etc.,  then  these  numbers  are  not  applicable;  others 
may  be  substituted. 

12.  This  is  an  example  (albeit  somewhat  detailed)  of  the  reasoning 
which  led  to  the  empirical  values  of  many  parameters  and 
functions  in  this  paper.  Refinement  of  such  estimates  based 
on  a  large  number  of  curves  is  necessary  prior  to  operational 
use  of  the  measures  generated  in  this  paper. 

13.  This,  of  course,  is  frequency  dependent. 

14.  If  the  two  curves  are  at  different  frequencies,  there  is 
considerable  latitude  in  the  choice  of  the  attenuation 
coefficient  a  for  this  curve.  It  appears  that  the  geometric 
mean  of  the  two  individual  attenuation  coefficients  is  a 
reasonable  value. 

15.  These  fluctuations  may  be  due  to  actual  phenomena,  model 
characteristics,  or  measurement  error. 

16.  In  this  step  we  know  that 

k2  <  01(il1)<  k2 

so  that  there  is  some  variation  between  T^  and  J^.  Consider 
three  possible  forms  for  T^(r)  whose  perturbation  function  P^(r) 
appear  as  the  three  curves  in  Figure  21.  Assume  is 

the  same  for  each  of  the  three  T^ ,  so  that  the  three  curves 
f^fr)  have  the  same  geometric  second  weighted  moments. 

^The  geometric  second  moment  of  P^(r)  with  weight  f(r)  is 
®l^'l^  *)  In  F^Sure  ^la  data  shows  fluctuations  of  period 
less  than  p,  such  as  could  be  produced  from  a  coherent 
summation  model.  In  this  case,  will  be  the  variation  from 
'top  to  bottom'  of  a  cycle  and  will  be  rather  small  compared 
to  p.  In  Figure  21b,  we  still  have  significant  variation 


with  a  distance  p,  but  will  now  be  close  to  p.  This 
behavior  (with  period  of  a  few  p)  represents  actual  conver¬ 
gence  zones,  where  one  can  drop  from  caustic  to  shadow  within 
distance  p.  In  Figure  21c,  we  see  behavior  of  a  period  much 
greater  than  p,  and  far  too  large  to  be  convergence  tones. 

This  sort  of  behavior  could  occur  when  a  propagation  loss 
curve  is  almost  monotonic,  but  is  essentially  20  log  r  for  the 
first  several  miles,  then  10  logr,  then  changes  due  to  the 
crossing  of  a  front,  or  a  rise  in  the  ocean  bottom.  A  single 
logarithmic  component  cannot  match  the  three  different 
types  of  behavior,  so  a  slowly-varying  function  is  produced 
for  p^.  In  this  case,  will  be  rather  small  in  comparison 
with  ✓e^T  •  This  is  the  rationale  behind  the  choice  of 
tests  in  the  second  step  of  the  convergence  zone  determination. 


To  summarize  this  three  step  test,  the  logic  is  given  in 
Figure  22. 


Particularly  those  with  coherent  phase  summation. 

This  is  equivalent  to  smoothing  the  transmission  loss  function. 

It  is  inadequate  to  simply  search  for  the  maximum,  for  that 
could  occur  at  range  6,  such  as  in  Figure  11. 

This  can  occur  when  the  second  derivative  changes  so  as  to 
indicate  a  perturbation  from  the  logarithmic  component.  In 
such  a  case  $k(R)  will  have  a  local  maximum,  but  T^(R)  will  not 

The  restriction  to  positive  X  should  be  obvious:  the  summation 
diverges  when  X  =  0 . 


As  before,  the  use  of  f(r)  as  a  weighting  function  is  to 
decrease  the  sensitivity  of  n  to  rv. 


c)  Slowly  varying  function,  indicating  no  convergence  zones 


SI 


Although  Q,  (r)  is  given  by  an  infinite  series,  there  are  at 


most 


e+r,  -R. 
k  1 


terms  which  are  not  identically  zero  on  the  interval 
[6,rk].  Thus,  the  evaluation  of  may  be  done  using  a 

finite  number  of  terms. 
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The  function  is  C  . 

If  one  formally  expresses  Q(0)  the  indeterminate  form  °°/°° 
leads  one  to  believe  that  ft(0)  should  have  the  value  Qk,°°. 

However,  this  is  not  the  limit  of  as  A^  decreases  to  0. 

Since  the  domain  of  A^  is  the  open  interval  (0,°°)  a  maximum 
need  not  exist  for  Q(Ak).  But  this  function  is  constant  for 
large  A  and  has  a  limit  for  small  A,  hence  the  supremum, 

which  always  exists,  is  either  the  maximum  (if  it  exists)  or 

the  limit,  i.e.  / 


sup  fi(Ak)  = 


max  fi(Ak)  if  it  exists 

lim  fi(A,  )  if  max  does  not  exist 

lv° 


Of  course,  whenever  this  limit  is  less  than  £ik,»  the  maximum 
exists. 

28.  If  the  maximum  peak  of  Pk  (as  opposed  to  T^)  occurs  at  R^, 
then  it  is  obvious  that 

C2(Ak)  <_  ^k,®  for  Ak>0 

and 

ft(Ak)  <  for  Ak<rk  +  e -  R1 

29.  The  desired  maximum  is 


Ak  -  R,  -  Rx 


>\1 

•  * 

'•A’l 


whereas  the  true  maximum  occurs  for  any 

>kirk+£  -Ri 

and  since  one  can  only  talk  about  a  second  convergence  zone 
range  if  3*2  is  less  than  r^,,  one  has 

Xk  -  rk  +  e  ”  R1  >  R2  +  £  ~  R1  >  R2  "  R1 

Therefore,  the  true  maximum  cannot  be  the  desired  maximum. 

Even  if  the  maximum  peak  of  is  not  at  (as  in  Figure 
13)  but  elsewhere,  trouble  can  develop.  In  Figure  14  the 
true  maximum  of  will  occur  at 

*k  =  R3  R1 

whereas  the  desired  local  maximum  is  at 
Xk  =  R2  "  R1 

I.e.  the  interval  should  not  contain  both  R^  and  R3  in  Figure  14 

Which,  of  course,  is  a  geometric  attribute  and  is  usually 
described  qualitatively. 

If  6/r^2  is  small  enough  (e.g.,  less  than  1/8)  then  if  the 
curves  are  close  their  coefficients  should  be  close  as  well, 
provided  that  the  attenuation  coefficient  a  is  the  same  for 
both  If  S/r^  lar9e  (e.g.  1/2)  or  the  two  curves  T^ 

are  for  different  frequencies  (so  that  the  attenuation 
coefficients  differ)  then  it  is  quite  common  to  have  the 
curves  relatively  close  over  [6,r12l»  but  have  the  coefficients 
bk,bk  vastly  different. 

There  is  always  the  possibility  that  one  curve  has  measure 

2 

just  below  the  threshold  k  ,  the  other  has  measure  just  above 
the  threshold,  and  the  two  curves  are  quite  close  but  the 


indicator  Tg  is  1.  This  is  a  problem  of  binary  decisions  in 
general.  The  decision  as  to  the  existance  convergence  zones 
(cf.  Fig.  22)  is  a  binary  choice,  even  though  transmission 
loss  curves  may  vary  continuously  from  a  no  convergence  zone 
situation  to  existance  of  convergence  zones. 

35.  One  may  also  devise  measures  based  on  Xi  <  such  as 

sC 

Xi  -  x2 

or 

X1  J  ®2(*2) 

X2  V  Q1(S,1) 

provided  x2  not  zero.  But  at  this  time  it  appears  that 
the  measures  and  will  be  adequate. 

36.  See  note  2. 

37.  Though  they  usually  will  have  a  peak  every  or  so. 

38.  See  note  5. 

39.  The  first  minor  is  obtained  by  multiplying  (-l)1+-,  by  the 

X.U.  i 

determinant  obtained  by  deleting  the  i  row  and  j '*  column 
of  the  determinant  A(k). 

40.  These  second  minors  are  obtained  by  multiplying  (-l)1+]+m+n 

i.  u 

by  the  determinant  obtained  by  deleting  the  i  and  j  rows 
and  mth  and  nth  columns  of  A(k). 

41.  The  reason  for  ignoring  the  reflected  paths  is  the  simplification 
of  the  analysis.  It  is  possible  to  do  the  analysis  with 
surface  and  bottom  reflected  rays,  but  then  the  field  would 

not  be  represented  as  easily.  The  family  of  ray  paths  would 
consist  not  simply  of  the  one  set  of  analytic  functions,  but 
would  be  a  countable  set  of  piecewise  analytic  functions. 


69 


«  •  •  -  •  .  • 


•  -  ..tV.  «*  - 


r.  ■'  >.  -r.  \  . 


*.\V 


42.  In  this  case 

_  9  Z  _  3  Z  o 

;  3  6.  3  e  d  0. 

1  o  l 

43.  If  it  were  zero  the  raypath  would  not  pass  through  the  source 

since  z,  4=  z  . 

1  o 

44.  As  a  result  of  that,  one  may  easily  obtain  branches  which  are 
false  roots.  kore  serious  for  numerical  work,  however  is  the 
problem  of  identification  of  branches,  since  a  single  value  of 
€  will  produce  eight  pairs  (x,y)  of  points  (not  all  of  which 
need  be  real) . 
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